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The early stage of a heavy ion cohision is governed by local 
non-equilibrium momentum distributions which have been approx- 
^ ■ imated by colliding nuclear matter configurations, i.e. by two 

CN , Lorentz elongated Fermi ellipsoids. This approach has been ex- 

pi2 I tended from the previous assumption of symmetric systems to 

\^ ' asymmetric 2-Fermi sphere configurations, i.e. to different densities. 

Cp . This provides a smoother transition from the limiting situation of 

two interpenetrating currents to an equilibrated system. The model 
is applied to the dynamical situations of heavy ion collisions at in- 
termediate energies within the framework of relativistic transport 
(RBUU) calculations. We find that the extended colliding nuclear 
^ ' matter approach is more appropriate to describe collective reaction 

^ . dynamics in terms of flow observables, in particular, for the elliptic 

flow at low energies. 
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1 Introduction 



Heavy ion collisions from Fermi energies up to relativistic energies of 1 — 2 
AGeV provide the unique possibility to explore the nuclear matter equation 
of state (EOS) under extreme conditions of density and temperature in the 
laboratory [1]. In a hydrodynamical picture the time evolution of such a re- 
action can be understood in terms of a pressure gradient which builds up in 
the compressed zone and drives the dynamics. Naively, one therefore expected 
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to obtain a direct access to the nuclear EOS in the measurement of collective 
flow observables. However, in fact, the situation is more complex: The system 
does not behave like an ideal fluid but binary nucleon-nucleon collisions lead 
to a viscous behavior. Moreover, calculations show that over most of the reac- 
tion time, in particular during the compressional phase, the system is out of 
/oca/ equilibrium, see e.g. [2]. Also experimental evidence for non-equilibrium 
in terms of incomplete stopping even in central reactions has recently been re- 
ported [3]. Thus the hydrodynamical limit is not reached in relativistic heavy 
ion reactions except in the final stages of the collision. 

On the other hand, microscopic transport models have proved to be an ade- 
quate tool for the description of the non-equilibrium reaction dynamics at in- 
termediate energies [1]. The physical input of such semi-classical models based 
on Boltzmann type equations are the nuclear mean field U and the nucleon- 
nucleon (NN) cross section a. Both are derived from the effective two-body 
interaction in the medium, i.e. the in-medium G-matrix; U ~ {^Gp), a ~ 
{'isGp), respectively da/dQ ~ ICp- However, in most practical applications 
phenomenological mean fields are used. Adjusting the known bulk properties 
of nuclear matter around saturation one has to rely on extrapolations to supra- 
normal densities. One can then try to constrain the models with the help of 
heavy ion reactions [4,5]. 

As mentioned above, in heavy ion collisions the system is far away from lo- 
cal equilibrium even during the high density phase which mainly governs the 
dynamics of the process. In a fully consistent treatment non-equilibrium ef- 
fects should be considered self consistently within the framework of micro- 
scopic many-body theory and a dynamical description of the reaction, e.g. 
by Boltzmann type transport equation [6]. An exact solution of the problem 
would require a consistent treatment of the transport equation and the mi- 
croscopic structure equations, e.g. within the framework of Dirac-Brueckner 
(DB) theory, for general non-equilibrium situations. Thus one would have to 
determine the relativistic in-medium G-matrix not only for an equilibrated 
Fermi sphere, but for arbritrary momentum distributions, in principle also at 
finite temperatures. This is presently not possible without introducing further 
approximations. The simplest approximation is the Local Density Approxi- 
mation (LDA), which assumes a local spherical, sharp Fermi sphere, and thus 
essentially neglects non-equilibrium effects in the in-medium interactions. As 
a step towards non-equilibrium the Colliding Nuclear Matter (CNM) model 
[7] was introduced (which was also called the Local phase space Configuration 
Approximation, LCA), which parametrizes the local momentum distribution 
by two Fermi spheres with a finite relative velocity. We have shown previ- 
ously that the non-equilibrium effects included in the CNM model modify the 
underlying EOS, in particular, the dynamical EOS effectively seen during a 
heavy ion collision is softened compared to ground state matter at equivalent 
densities [8]. This affects the collective dynamics, e.g. in terms of observable 



flow signals [9,10]. From these studies it was concluded that this type of non- 
equilibrium features should be included at the mean-field level in transport 
analysis of heavy ion reactions in order to allow more reliable conclusions on 
the nuclear matter EOS. 

The Colliding Nuclear Matter (CNM) model has been extensively studied in 
[7] and has been applied to intermediate energy heavy ion collisions in refs. 
[8-10]. There it was shown that the consideration of non-equilibrium effects 
within the CNM approach essentially changes the description of collective 
dynamics relative to a treatment in the local density approximation (LDA). 
Differences of the two approximations appeared in collective flow observables 
and were found to be of the same magnitude as those which arise due to a 
different, i.e. soft or a hard, EOS. 

In this work we extend our previous studies by including more details in the ap- 
proximations involved in the description of non-equilibrium situations within 
the CNM approach. In previous applications a symmetric conflguration was 
assumed, i.e. two nuclear matter streams of equal (averaged) density. This 
should not be very realistic, particularly in peripheral collisions. For this pur- 
pose the Asymmetric Colliding Nuclear Matter (ACNM) model is introduced 
here, which accounts also for a possible asymmetry of the Fermi momenta of 
the two Fermi fluids. This situation arises in those regions where the tail of one 
nucleus penetrates into the interior of the other (see e.g. Fig. 1). In previous 
studies based on the CNM approximation we mainly considered symmetric 
colliding systems. A locally asymmetric two-stream situation is also more fre- 
quently encountered in collisions of different mass systems. Such reactions 
will experimentally be studied more often in the future. Thus the extended 
approximation is a further step into the direction of ref. [6], i.e. the attempt to 
describe the complex dynamical situation in heavy ion reactions as precisely as 
possible. The extension to ACNM is found to be particularly important at low 
and intermediate energies when peripheral heavy ion collisions are considered. 



2 Approximations in colliding nuclear matter 



As discussed in the introduction the treatment of non-equilibrium effects in 
dynamical transport descriptions cannot be done without any approximations 
due to its high complexity. In the CNM approach [7] ground state DB re- 
sults are extrapolated for idealized 2-Fermi spheres or, covariantly, for 2-Fermi 
ellipsoids, at zero temperature. This approximation constructs the 2-Fermi 
sphere conflguration in a covariantly consistent way with respect to the effec- 
tive masses, but neglects the blocking in the intermediate propagator due to 
the second current. It was shown in [7] that this approximation is reliable due 
to a moderate dependence of the DB self energy around the Fermi momentum. 



In applications to heavy ion reactions the relevant approximations are (a) 
the zero-temperature limit and (b) the restriction to symmetric CNM distri- 
butions, i.e. same Fermi momenta of the two counter-streaming currents. To 
demonstrate these effects we show in Fig. 1 typical local momentum distri- 
butions as they are obtained from transport descriptions of central Au+Au 
reaction at 0.6 AGeV beam energy at different locations in coordinate space 
(in the center of mass, and in beam direction) at the time just before the com- 
pression phase. In the left column we show the local momentum distributions 
as they are obtained from relativistic BUU calculations. Details of the RBUU 
model used here are given in section 4. In the right column are asymmetric 
fits to the momentum distribution, i.e. fits with two Fermi distributions (in- 
cluding Pauli corrections) of different density and temperature, for details see 
[2]. It is seen from these fits that the momentum space consists in general of 
two rather different Fermi distributions outside of the central cell, which are 
furthermore characterized by different local temperatures of the subsystems. 



How to go beyond the zero-temperature approximation is difficult to study. 
The DB approach [11,12] has been mainly applied for ground state nuclear 
matter at zero temperature. The few existing finite temperature calculations 
indicate a moderate T-dependence of the mean field at temperatures of a sev- 
eral MeV [11]. Thus, one can expect that the approximation should be reliable 
for small thermal excitations with moderate momentum tails. In heavy ion col- 
lisions high thermal excitations with local temperatures of T ~ 20 — 40 MeV 
in the central region during the compression phase can be reached at inter- 
mediate energies, (see Fig. 1), which makes this approximation less obvious. 
However, we believe that temperature effects become important in a higher 
order correction and, in addition, complicate the calculation significantly. The 
determination of a local temperature in heavy ion collisions, which is necessary 
for such a treatment, requires a fit procedure at each space-time point. Thus, 
we do not consider a possible temperature dependence of the mean field here. 



The assumption of a symmetric CNM approximation should be reliable for 
central collisions and high beam energies. However, as seen in Fig. 1 this 
is not the case especially if one goes away from the central cell. We expect 
the asymmetry of the CNM configuration to be particularly important at 
low energies. At high energies, on the other side, faster time scales and a 
higher transparency of the system - for peripheral reactions - leads rather to a 
separation between 2- and 1-Fermi sphere configurations. Thus, we extend the 
CNM model to asymmetric 2-Fermi ellipsoid situations (called as Asymmetric 
Colhding Nuclear Matter (ACNM) in the following). 



Asymmetric colliding nuclear matter in the relativistic DB ap- 
proach 



The treatment of non-equilibrium effects within the spirit of a CNM approx- 
imation has been extensively investigated in ref. [7]. Here we extend this ap- 
proach to Asymmetric Colliding Nuclear Matter (ACNM) by introducing an 
asymmetry in the Fermi momenta of the two counter-streaming systems. We 
describe the idea of the formalism here; more detail can be found in ref. [7]. 

The momentum distributions are described in terms of a superposition of two 
counter-streaming nuclear matter currents, i.e. two boosted Fermi ellipsoids 
at zero temperature. This configuration is schematically shown in Fig. 2. An- 
alytically we write it as 

/l2(k) = /i(k, kp,) + UK kp,) + (5/(k) 

= 0(^t - KO + ^{^^l - ^X) + '^/(k*) . (1) 

O is the step function, kp^ (i=l,2) are the Fermi momenta and u'( = (71, Ui7i) 
are the streaming four velocities of the two boosted nuclear matter currents. 
The chemical potentials in the moving system are then given as 



fii{k) = fx*{k) + S„< = ^k%^ + m*2 + ^^u^ 

in terms of the effective chemical potentials fi*{k) in the rest system. The last 
term in eq. (1) Sf{k) guarantees the conservation of the Pauli principle in the 
overlap region between the two nuclear matter currents, and it is given by 
(5/(k) = — 0i(k)G2(k). However, to guarantee baryon number conservation 
one has to restore the total baryon density. This can be done in a covariant 
manner by redefining the Fermi momenta [7]. 

The ACNM configuration is characterized by three invariant parameters, the 
invariant densities of the two Fermi ellipsoids and the relative velocity 



Pi = \/Jif^Jt > ^rcl 



Vi - V2 



V1V2 



(2) 



An alternate set of configuration parameters, denoted collectively by X; ell's 
X = {vrei, Ptot, Ps} where ptot = Pi + P2 is the total invariant density, and 

Pi - P2 .„^ 

ps = (3) 

Ptot 

the density asymmetry parameter, ps = corresponds to the special cases 
of symmetric CNM, p^ = 1 to a single Fermi-sphere at rest with invariant 



density pi = ptot and a single nucleon with relative velocity frci- Thus, the 
ACNM model naturally includes a smooth transition to the local density ap- 
proximation (LDA) of equilibrated nuclear matter for both cases f rei *— ^ and 
p^ I— ^ while symmetric colliding nuclear matter contains only the first limit. 
All the expressions for the scalar and baryon densities can be derived in terms 
of these invariants. 

The relativistic self-energy in CNM as well as in ACNM has the same Lorentz 
structure as in ground state matter, i.e. it contains a scalar part entering 
into the effective mass m* and a vector part entering into the kinetic four- 
momentum k* [7,8] 



j:{k) = j:s{k)-j,^^{k) . (4) 

The construction of ACNM is done in a manifestly covariant way and it can 
be defined in any reference frame. It is, however, most appropriate to work 
locally in that frame where the total baryon current vanishes ji2 = ji +32 = 0. 
We will denote it as current-zero frame RSu in the following and we will 
work always in RSu- The space- like part of the vector field vanishes in -RS'12 
by definition S12 = Si + Si = and k* = k. This simplifies the situation 
considerably, since canonical and kinetic momenta are not identical in general 
frames. 

As discussed above, we construct ACNM by extrapolating from calculations of 
ground state nuclear matter in the DB approach. The self-energies are thereby 
given as 

^m{k; kp) = C J d^q gm{\f\\) /(q; kp) %n{k, q) (5) 



in terms of the distribution function /(q; kp) of one Fermi sphere (C = -n^^ 
with K = 4 for nuclear matter). Here m = s, stands for scalar and (time 
component) vector self-energies. The T matrix amplitudes can be taken from 
DB calculations and Qm in (8) are appropriate weight factors corresponding 
to the Lorentz structure of these amplitudes [7,8]. The invariant T matrix 
amplitudes themselves are determined in equilibrated nuclear matter, i.e. for 
a single Fermi sphere. We further define effective coupling functions as 



In Relativistic Mean Field (RMF) theory these quantities correspond to glnlrv?^^ 
where Qm are the coupling constants of the m-meson to the nucleon and m^ 
its mass. 



The ACNM self-energies are given in a corresponding way by using the ACNM 
distribution function of eq. (1) for a configuration specified by the parameters 

X = {^rchPtot,^^} 



j:';^'\k;x)=C J d'qgU\q\)fi2{q;x)'rUk,q,x) 

= ^^i\k-x) + ^^Z\k-x) + 5^^{k-x) ■ (7) 

In the second hne it is written as a sum of contributions of two streaming 
Fermi spheres and a Pauh correction. Note, that the effective T matrix also 
depends on the configuration parameters, since the intermediate propagator 
does. In the first two terms the self energies for the two streaming nuclear 
matters are specified as (i=l,2) 



E«(A;;x)=c|rf3g^„(|q|)/i(q;;^)T„(A;,g,x) 



The second line contains the essential approximation of the approach, namely 
that the coupling function of ground state nuclear matter are used at momenta 
corresponding to the appropriate Lorentz boost multiplied by the density of 
the ACNM configuration. The Pauli correction can be evaluated in a similar 
way. Since the fields are only moderately momentum dependent below the 
Fermi surface, the coupling functions can be determined by evaluating at the 
Fermi momentum. For the special case of symmetric matter (CNM, kp^ = 
kp^ = kp) one obtains 

ST^mik; X) = C j (fq fi-rndql) Tm{k, q) 6f{q; x) = Tm{kp, k)5pm{,x) 

^Tmikp,kp)5pra{x) ■ (9) 

In this form all the self energies depend on the variables {k; x} = {k; frci, Ptot, Ps}- 
They are difficult to handle in a transport calculation, in particular because 
of the explicit momentum dependence. In analogy to the Hartree (or RMF) 
approximation we want to obtain momentum-independent self energies. These 
are obtained by averaging over the momentum of the ACNM configuration 

^(12), , _ Jd^kgm{H)i:g^^{k,x)fi2{Kx) _ r(i2). X n2).x ..n^ 
^m {X — Tyrj 71 IN r n ^ — ■*- m {X Pm {X A-i-UJ 

where analogously as in eq. (6) ACNM coupling functions are defined in the 
second equality. Details of the averaging procedure and the calculation of the 
coupling functions are given in the Appendix. 



Fig. 3 shows the density asymmetry dependence of the effective couphng func- 
tions r^Q (see Eq. (10)) at fixed relative velocity and total densities (solid 
circles). It is seen that the effective couplings rise with increasing asymmetry 
ps and approach the LDA limit. 

At high relative velocities and moderate densities, where Pauli effects are only 
of minor importance, a rising asymmetry parameter shifts the center-of-mass 
towards the current with the higher density and the mean field is dominated 
by the momentum distribution of this ellipsoid. Thus, in averaging over the 
explicit momentum dependence of the fields, mainly moderate relative mo- 
menta contribute to the ACNM self energy. On the other hand, the invariants 
Tsfi increase with decreasing momentum which leads finally to the observed ps 
dependence of the Fig. 3. This behavior is more pronounced at larger densities 
since there the variation of the underlying T-matrix amplitudes is stronger. 

The ACNM mean field depends on three parameters which complicates their 
application in transport calculations for heavy ion collisions. For practical use 
we apply therefore a simple parameterization of the form ([/(p, fj-ei, (5) stands 
for scalar and vector self energies) 

U^^^^ ^U{p,v,^,,ps)= p]U'^^\p) + {l- p])U'^''''{p,v,^{) . (11) 



This parameterization is also shown in Fig. 3 (solid lines). It is sufficiently 
accurate to be used in heavy ion collisions. 



4 Application to heavy ion collisions 



For the theoretical description of heavy ion collisions we use the relativistic 
Boltzmann-Uehling-Uhlenbeck (RBUU) transport equation [6,10]. The mean 
field will be implemented in three approximations as discussed in the previous 
section: The symmetric CNM model has been studied in detail previously 
[9,10], its extension in the ACNM approach based on the parameterization 
(11) is studied here for the first time. Results are also compared to the LDA 
where the mean field is only a function of the total density, i.e. taken at 
zero relative velocity Vre\ = 0. The density and momentum dependence of 
the mean field is taken from the DB model of Ref. [12]. The nuclear matter 
saturation properties obtained with the Bonn A potential are psat = 0.185 
fm~^ and E = —16.1 MeV. With a compression modulus K = 230 MeV 
one has a relatively soft EOS. In this context one should keep note that the 
consideration of the non-equilibrium effects in the CNM/ACNM approaches 
weakens the nucleon effective interaction leading to an even softer EOS than 
that given in the LDA. In all the simulations the collision integral is treated in 



the same way by using standard parameterizations for the elastic and inelastic 
total and differential cross sections [13] including A and A^* resonances. The 
resonances are propagated in the same mean field as the nucleons and they 
decay into one and two-pion final channels. The pions are propagated under 
the influence of the Coulomb interaction. However, they interact strongly with 
the hadronic environment due to absorption. 

We have analyzed Au+Au collisions at intermediate energies in terms of col- 
lective flow effects given by its in-plane and out-of-plane components. They 
can be characterized by the first and second order Fourier coefficients of the 
azimuthal distributions, N{(f)) ~ l+ficos(0) + 2f2cos(20) + - ■ -, and calculated 

as 



In general, f 1^2 depend, apart from the beam energy and centrality, on rapidity 
and transverse momentum. This has been studied extensively experimentally, 
e.g., in refs. [1,14], and theoretically using non-relativistic approaches for the 
mean field in [4,5,15-18], and relativistically in [9,10,19] and [20] where in the 
latter case it was claimed to see a softening of the nuclear EOS as nuclear 
density becomes large. 

We begin the flow analysis by discussing the rapidity and transverse momen- 
tum dependence of the in-plane flow vi in the Figs. 4,5 for semi-central and 
peripheral Au-|-Au reactions at 0.25 and 0.4 AGeV. For a realistic simulation 
of the experimental conditions the RBUU events were passed through a phase 
space fragment coalescence algorithm including a reaction plane determina- 
tion in the same way as the FOPI experiment [21,22]. The centrality selection 
was performed using the total multiplicity of charged particles, see e.g. [10]. 

From these figures it can be seen that the non-equilibrium effects accounted by 
the ACNM and CNM models influence the directed in-plane flow only mod- 
erately. The differences are relatively biggest at large transverse momenta p\ 
and in the high rapidity regions. This is expected, since high p\ particles are 
mainly emitted earlier when the matter is still anisotropic. Mostly, the in-plane 
flow with ACNM mean fields lies between the CNM and LDA results. This 
result is to be expected since ACNM interpolates the self-energies between 
the two extreme cases, CNM and LDA. The different flow pattern originate 
from the treatment of momentum anisotropies on the mean field level. The 
momentum dependence included in the CNM and ACNM approximation gen- 
erally weakens the repulsive vector self energy, in particular at high incident 
energies, which leads to an effectively softer EOS compared to the nuclear 
matter EOS (LDA case) [8]. This mechanism reduces the in-plane flow. 



The differences become more pronounced for the elhptic ffow shown in Figs. 6,7. 
The elhptic flow is generally considered as a suitable tool to study the momen- 
tum dependence of the nuclear mean field [14,4,25]. The elliptic flow is created 
during the high compression phase and governed by the pressure gradient due 
to the expansion of the initially compressed fireball as well as by shadowing 
effects of the spectator matter. Since the high density phase of the collision is 
still governed by local non-equilibrium [8] one expects to see clear differences 
between the CNM, ACNM and LDA models. The model dependencies are 
most pronounced in the transverse momentum dependence of the elliptic flow 
(Fig. 6). The energy dependence of the p\ integrated elliptic flow (Fig. 7) 
shows also clear signals from the highly anisotropic compression phase. The 
LDA yields a very strong elliptic flow and fails to describe the excitation func- 
tion above 0.25 AGeV. The CNM approximation, on the other hand, reduces 
the flow, particularly, at high energies and provides a good fit to the data 
there, but is too small at low energies. At low energies the ACNM model lies 
between the two limiting cases, LDA and CNM, which results in a rather ac- 
curate description of the excitation function over the complete energy range 
from 0.1 H- 4 AGeV. 

The differences between LDA, CNM, and ACNM approaches have the same 
origin as for the in-plane flow although the model dependencies are now more 
pronounced. Accounting for momentum anisotropies on the mean field level 
the repulsion of the model is weakened. Thus, due to the smaller pressure gra- 
dients the initially compressed system expands slower which effectively reduces 
the shadowing effect of spectator matter and the magnitude of the out-of-plane 
emission. The differences between CNM and ACNM can be understood from 
Fig. 3. We also find that the asymmetry dependence is similar at other rela- 
tive velocities and densities not shown here. The increase of the vector field 
with asymmetry parameter ps result in a stronger elliptic flow in ACNM as 
compared to CNM, but still smaller with respect to LDA, as expected. One 
obtains a smoother transition to LDA as the beam energy decreases. 



5 Summary and conclusions 



We studied heavy ion collisions with a relativistic transport model and fo- 
cused on non-equilibrium features of the phase space and their correlation to 
collective dynamics in terms of flow signals. The non-equilibrium effects were 
described by a Colliding Nuclear Matter model. We extended this model to two 
Fermi ellipsoids with an asymmetry in their densities, called as ACNM. The 
application of CNM and ACNM at intermediate energy heavy ion collisions 
was compared to a simple Local Density Approximation (LDA). 

We analyzed collective flow effects in and out of the reaction plane. It was 
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found that the directed in-plane flow is only moderately affected by asym- 
metry effects in the ACNM approach, except if one considers the transverse 
momentum dependence at high pf. For high energetic particles the ACNM 
model provided the best description of the data especially at the lowest bom- 
barding energy considered (0.25 AGeV). 

The elliptic flow excitation function is more sensitive to these different ap- 
proaches since this observable is built during the early non-equilibrium stage 
of the collision. We observed a reduction of the elliptic flow in the ACNM 
and CNM models compared to LDA, where the ACNM model gave the best 
description of the data at intermediate energies. The asymmetry of the conflg- 
uration was found to be particularly important at low energies and provided 
a smoother transition to the LDA. 

We interpreted the results by the fact that the repulsive vector fleld is essen- 
tially reduced in the CNM and ACNM cases which also decrease the magni- 
tude of flows as compared to LDA. This considerably improves the comparison 
with flow data. The deviations between CNM and experiment especially at 
low energies in the elliptic flow excitation function seems to be resolved when 
asymmetric CNM conflgurations are considered (ACNM) where a smoother 
transition to the LDA case is now observed. 

We conclude that non-equilibrium effects are important describing the collec- 
tive dynamics of intermediate energy heavy ion collisions. A fully consistent 
treatment of heavy ion dynamics and nuclear structure equations within the 
framework of non-equilibrium transport theory is presently not possible. The 
CNM/ACNM models provide an approximative, but solvable tool to incor- 
porate local non-equilibrium features more consistently at the level of the 
effective interaction. For deflnitive conclusions on the nuclear matter EOS it 
is important to study such approximations in complex systems. Our results 
can be considered as a flrst step in this direction. The next step would be to 
consider flnite temperature effects, as well as the explicit momentum depen- 
dence of the flelds. 
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Fig. 1. Left: Local momentum distributions from RBUU at positions 
X = {x,y,z) = (0,0,0) (top), X = (0,0,2) (middle) and x = (0,0,4) (bottom) 
(in units of [fm]) at time t = 10 fm/c. The figures on the right give the corre- 
sponding ACNM momentum distributions at finite temperatures Ti^2 with values 
of ri,2 = 29 MeV (top), Ti = 40 and Ts = 8 MeV (middle) and Ti = 25 und 
T2 = 1 (bottom). The considered reaction is a central (6 = fm) Au+Au reaction 
at Eheam = 600 AMeV beam energy. 



12 



NMl 




Vrel 



Fig. 2. Schematic representation of the ACNM configuration given by two co- 
variant Fermi ehipsoidal momentum distributions /i, (i=l,2). These are charac- 
terized by Fermi momenta (6i = kp^) which are elongated along the boost direction 
(ai = 7ififcFi). 
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Fig. 3. Effective coupling functions F^ q (left and right figures respectively) as 
function of the asymmetry parameter ps at fixed total densities of ptot = 1 [Psat] 
(top) and ptot = 2 [psat] (bottom) and relative velocity v^ei = 0.5 [c]. The circles are 
results of ACNM calculations and the solid lines indicate the parameterization of 
Eq. (11). 
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Fig. 4. In-plane collective flow v\ for charged (Z = 1) particles in Au+Au reactions 
at 0.25 AGeV incident energy as function of the normalized transverse momentum 
'PI at different rapidity intervals. The centrality intervals M4 and M3 corresponds 
to impact parameter ranges of (3 — 6 fm) and (6 — 7.5 fm), respectively. RBUU cal- 
culations treat the mean field in the following approximations: local density approx- 
imation (LDA) where no configuration dependence is included (dot-dashed lines); 
colliding nuclear matter (CNM) approximation (dashed lines); asymmetric colliding 
nuclear matter (ACNM) approximation (solid lines). The FOPI data (filled circles) 
are taken from [22]. 
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A Appendix 



The ACNM configuration (1) given by two counterstreaming nuclear matter 
currents is characterized by three parameteres, namely the Fermi momenta 
kFi2y respectively the rest densities of the currents, and their relative velocity 
frei- Expression (1) contains a Pauli correction term 6f = — 0iO2 which re- 
stores the Pauli principle in case that the relative velocity of the currents is 
small and the original ellipsoids overlap. As discussed in detail in [7] the Pauli 
correction can be performed in a covariant way. Corresponding expressions for 
the Pauli correction terms can be found in [7]. 

In [7] the mean field, respectively the counterstreaming self-energy components 
have been derived for symmetric nuclear matter configurations {kpi = kp^). 
Here we present the extension to the general case of asymmetric configurations 
{kp^ 7^ ^^2)- In order to obtain mean fields in a Hartree form which depend 
only on the ACNM parameters {kp^^^v^^h Ps] the self energy contributions of 
the two currents S^\ (i = 1, 2, (5, m = s, 0) have to be averaged over the ACNM 
configuration in momentum space. This includes also corrections arising from 
the conservation of the Pauli principle. 

Since the underlying DB self-energies show only a moderate momentum de- 
pendence below the Fermi surface {k < kp) [12] the integrations over the 
comoving current and also over the overlapp region are numerically easy to 
carry out due to moderate relative momenta. The difficulty appears for the re- 
maining integrations of the self-energies over the second Fermi ellipsoid (i 7^ j) 



V(i) f. 



It is convinient to perform the integrations of the self energies S^^ in their cor- 
responding rest frames RS\. In this frame they depend only on the modulus of 
the 3-momentum which reduces the integrations (in polar coordinates) to two 
dimensions. In particular, one has to perform integrations over spherical cells 
which include the integration area over the second ellipsoid as schematically 
depicted in Fig. A.l. 

Let us for example consider the integration of S^^)|/?5^ over the other ellipsoid 
5*2, which has to be performed in the rest frame RSi (Fig. A.l with i=l, j=2): 



\/k'^ + m* 



d?kTsi{kp,,k)h — = m*C J edk j dn{k)Tsi{kp„k)f2 



n{k) 




Fig. A.l. Schematic representation of the averaging of the self energies over the 
ACNM distribution. The integration of Sm|/j5. has to be performed in the rest 
frame RSi by integrating over the other elhpsoid in polar coordinates. 

The integration of the spatial angel cos 6 = -|- yields 



k{:'l 



I dQ{k) = -2n I = 27r (-1 - -^ 



n{k) 



ki'^/k 



The intersection point k^f"^ of the sphere with radius k^ with the second ellip- 
soid 2 is thereby given as [7] 



ki^^ 






(A.l) 



Using eq. (A.l) with Vi = 0, f2 = — frd, Q'l = 0, 71 = 1 und kp-^ = k one 
obtains 



Kz \RSi 



Ef2 - 7(^rci) Vfc^ + m*'^ 



In summary the integration over the polar angles leads to 

,---,-, ,, -VWTm^)] in RSi 




VW+m^)] in RS2 



In general, there exist different cases depending on the position of the center 
of the ellipsoid Sj, i.e. it can lie inside or outside the sphere with radius k. 
This leads to different cases for the integration / dfl: 
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S] outside of RSi ■ 



Sj inside RSi 



{ \qj\ -ttj <k< \qj\ + as J dn{k) = 27r(l + A;W/A;) 
I otherwise j dVL{k) = 

()<k<aj-\qj\ JdQ{k)=ATT 

aj - \qj\ <k <aj + \qj\ J dn{k) = 27r(l + ki'^k) 



(A.2) 



k > ttj + \qj\ 



Jdn{k) =0 



The cases (A.2) can be summarized. For the scalar part of the self energy T,^^ 
the integration over the second ellipsoid Sj in the rest frame RSi reads 



Umax 

m*C I k'dk 

k 



2 ,, ^Akp^k) 



\/k'^ + m 



*2 



dViik) 



n{k) 



max{0,aj — \qj\} 



27rm* ■ 



Tsi{kF,,k) ^,^2 



y/W^+rn*^ 



2k'dk + 



Tsi{kF,,k) 



l«j-kjl 



V^^ + 1Ti 



*2 



e + 



k (Et 



Vrel V 7 



— ykP^+mF^ 



dk 



The total scalar part of the ACNM self energy averaged over the ACNM 
configuration reads finally 



rnax[0,'y{kp^-Vr^lFp^)\ 



27rm*C 



Ts2{kF2,k) ^,2 



■yi^kF-^+VrelEp-^) 



2Txm*C 



^Ikp-i^-VrelEp-^ 

. max[0,-y{kp^-v^^iFp^)] 



Ts2ikF2,k) I 



m 



*2 



2k^dk + 



k (E 



k-^^—it^- v/pT 



m 



*2 



Vrel \ 7 



Tsl{kF^,k) ^2 



yW^+wF^ 



l[kF2+'"relEF2 J 



Tsi{kFi,k) 



2k^dk + 



k (E 



^2 ^ j;_ -^ _ ^/^^ 



'-l\kp^-VreAEpA 



\Jk'^ + m*2 \ t;re/ V 7 



m 



*2 



2<^Ps(x) Y^ ^sjikF^)ps\x)+'^s{min{kF„kF2))5p1{x) 



UA; 



pi^Hx) + 



dk 



Pl'^ix) + 



(A.3) 
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The scalar densities depend on the ACNM parameters x due to the config- 
uration dependence of the effective mass m* = m*{x)- The determination of 
the vector part of the mean field can be performed in a similar way. First one 
defines it covariantly by 



= < S(12)^(A;; X)fl2{k; X)^ > Jl2./jl2aj'^2 

^TZjl2./jl2aJ^2 , (A.4) 

with T^^ =< Tj^^"^^^ fi2k'^ / E >. Eq. (A.4) is valid in any reference frame, hence 
also in the special frame RS12 where the total baryon current vanishes: 



Toi U^Sf + Sff,) + ro2 i&f + 6ff,) + 

r05Sj^Sf]jl2u/jl2aJ?2 ■ (A.5) 

Here Foj contains the integration over the comoving ellipsoid RSi, and Fqj 
denotes the corresponding integration over the second ellipsoid j. These are 
exactly the same integrations as discussed above. 

From the averaged self energies (A. 3, A. 5) one can finally derive effective cou- 
pling functions for the ACNM configuration according to Eq. (6) with m = s,0 

^m (X) = (12) • (A.6) 

Pm [X) 
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